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Excitation Mechanisms of the Oxygen 5577 Emission 

in the Upper Atmosphere 

E. Tandberg-Hanssen and F. E. Roach 

(June 18, 1959) 

Possible excitation mechanisms for the green 5577 emission are considered in the light 
of recent data on the dynamics of the upper atmosphere. Photochemical reactions as 
affected by mass motions as well as excitation directly due to the mass motions are analyzed. 
It is concluded that either or both mechanisms could probably account for the observed 
emission. 



1. Introduction 

It is now well established (hat the hulk of the 
X5577 airglow is generated in a fairly thin layer of 
the atmosphere at a height of about 100 km (Hepp- 
ner, Stolarik, and Meredith [l]; 2 Roach, Megill, 
Rees, and Marovich [2]; and Tousey [3]). Further- 
more, it is known that the airflow 3 is not emitted 
in a uniform static layer, but that the emission is 
concentrated in cells that move relative to the earth's 
surface (Roach, Tandberg-Hanssen, and Mem II 
[4,5]. 

A general theory of the 5577 nighl airglow which 
includes both the microscopic (excitation) mecha- 
nism and the macroscopic (dynamical cells) features 
is required, and in this paper two different possi- 
bilities are examined. 

The principal observational facts any theory will 
have to explain are summarized: 

(1) The average night intensity of the 5577 emis- 
sion corresponds to approximately 2xl0 8 transi- 
tions/cm 2 column sec. For a night of roughly 10 hr 
duration, 2 X 10 8 X 3.6 X 10 4 = 7 X 10 12 * transi- 
tions/cm 2 column are consequently needed. 

(2) The airglow intensity varies during one night 
sometimes by as much as a factor 5, and the maxi- 
mum comes generally not in the beginning of the 
night (Roach, Tandberg-Hanssen, and Megill [4]). 

(3) The emission is not uniform over the sky, but 
appears in "cells" having horizontal motions of the 
order of 100 m/sec. 

The energy necessary to excite the *S level of the 
oxygen atom (from which the atom may jump to 
the l T> state and emit the green line) is 4.2 ev (32,000° 
K). The temperature in the 100-km region is less 
than 300° K which corresponds to 0.04 ev, so it is 
impossible that any significant excitation can be 
brought about by the thermal motions of the par- 



i High Altitude Observatory, University of Colorado and National Bureau 
of Standards, Boulder Laboratories. 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 The following is restricted uniquely to the green emission, X5577. 



tides. A static theory must therefore draw on 
exothermic photochemical reactions which are con- 
sidered in section 2. Similarly, a dynamic theory 
must invoke environmental effects whereby energy 
is given perferentially to selected particles, as 
discussed in section 4. 

2. Photochemical Reactions 

Above 100 km most atmospheric oxygen is in the 
atomic form, and according to Scaton [6], the most 
important recombination mechanism will be of the 
3-body type, as proposed by Chapman [7]. This 
conclusion is adopted here and the Chapman reaction 
is considered as the most favorable photochemical 
excitation mechanism for the green oxygen line. 
The reaction may be written 

O+O+O-^OOSHOo. (1) 

The reaction coefficient k for reaction (1) may 
probably be written (Bates [8]) 

k= 5 X 10~ 34 Ti cm 6 /sec. (2) 

There is, however, no direct measurement of k and 
the dependence T% is purely conventional. The 
yield of oxvgen atoms excited to the J S state is given 

by 

F( 1 S) =^(0)] 3 cm" 3 sec" 1 (3) 

where p is an "efficiency factor/' giving the fraction 
of 3-body collisions (1) that lead to excitation of the 
*S state, and n(0) is the number density of atomic 
oxygen. 

Since the distribution of excited atoms and mole- 
cules over the different energy levels is not known, 
the value of p is not known. In these numerical 
calculations, the simplifying assumption is made that 
p=l/n where n is the number of different types of 
reactions (1). Taking p=0.2 may be a large 
efficiency factor, but it is at present probably the 
best. 
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Because of the 3-body collision, the airglow 
emission will be critically height dependent if the 
distribution of atomic oxygen is at all sharply peaked 
at any height. If r c is the lifetime of the *S level for 
collisional deactivation and r r is the radiative life- 
time, then the number of photons emitted at any 
level by metastable oxygen atoms by (1) is propor- 
tional (Bates and Massey [9]) to 



NO)] 3 ? 

Tc + Tr 



(4) 



Detailed calculations show (Nicolet [10]) that n(0) 
is peaked near 100 km. According to (2) one should 
hence expect a pronounced maximum in the green 
emission near 100 km, and this is in agreement with 
observations. The number of atoms formed per 
square centimeter column per second in the *S state 
is then given by 



Jo Jo 



(5a) 



To take into account the actual distribution of n(0) 
with height may be quite complicated. On the other 
hand, it is known that it is sharply peaked at £=100 
km, so most of the contribution from the integral in 
eq (5) com.es from that height. The distribution can 
hence probably be well approximated with a cosine 
function or an exponential. An exponential model 
with a peak value of n(0), taken from Nicolet 's [10] 
work, is used here. Putting 



n(0)=n max (0)e 



-ZJH 



where H is the scale height at 100 km, we find 



N( 1 S)=pk[n D 



M? fo 



z,H dz 



= l/3K[n max (0)fH=l/3Y Q [(^)H;K=pk. (5b) 

Now, according to Garstang [11], the radiative 
lifetime r r for the *S state is 0.74 sec, and it is con- 
ceivable that collisional deactivation could take 
place before the NQS) atoms jumped to a lower 
state and emitted the green line. Calling the proba- 
bility of deactivation d 3 , the number of transitions 
1 S— >*D is given by 

N(5577) = l/3Y ( 1 S)H(l-d 3 ) 

transitions/cm 2 column sec. (6) 

The different quantities entering eq (6) are dis- 
cussed in order to see what numerical values for the 
green intensity are to be expected from the Chapman 
reaction. 

(a) Reaction coefficient, k. 

There is no evidence of the effectiveness of the 
Chapman reaction (see Bates [8]) so the best is to 
suppose that it is of the same order of magnitude as 



for other reactions of this type. 

For r=200° K, eq (2) gives £~7X10- 33 

(b) Density of atomic oxygen, n (0). 

The best estimates (see Nicolet [10]) indicate that 
% ax (0) = 1.5xl0 12 cm -3 between 90 and 105 km, and 
the density drops off sharply below and above this 
region. 

(c) Scale height, H. 

The hydrostatic scale height H=kT/mg at this 
height in the atmosphere is H=6.5 km (U.S. Exten- 
sion to the ICAO Standard Atmosphere [12]). 

(d) Probability of deactivation, r/ 3 . 

The problem of deactivation of the 2 S level has 
mainly been considered in connection with the in- 
tensity ratio of the red to the green emission under 
auroral conditions. In view of the comparatively 
short radiative lifetime of the *S level, it has generally 
been assumed that collisional deactivation of this 
level may be neglected, i.e., d 3 =0 (see Seaton [13]). 

Whereas deactivation by electron collisions may 
be of importance for bright auroral displays (Seaton 
[13]), the electron concentrations required seem to be 
too high for the process to be of importance under 
ordinary airglow conditions. The rate of recom- 
bination (0+0^0 2 ) is so slow that there is a negligible 
decrease in the atomic concentration during the 
night (Bates and Massey [9]). 

In view of the uncertainties in the data, a deacti- 
vation of the *S level cannot be ruled out. However, 
it is felt here that there is no obvious process by 
which the *S level will be collisionally depopulated. 
The factor (1— rf 3 ) is therefore probably close to 
unity, and this value will be adopted here. 

Numerical values for the quantities in (6) will now 
be put in to see what intensities are to be expected for 
the green airglow according to the 3-body collision 
theory. Figure 1 shows the intensities Q in ray- 
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Figure 1. Predicted intensity of [OI] 5577 from eq (5b). 
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leighs 4 as a function of K for different values of n (0), 
boxing in the probable range of these parameters. 
It is seen that already the minimum probable values 
give a very bright airglow. Either the yield of 
excited oxygen atoms, (3), has been overestimated, 
or the importance of collisional deactivation has 
been underestimated. If neither is the case, it is safe 
to state that within the accuracy of these data, the 
Chapman reaction can not only explain the observed 
airglow intensities, ,but it actually predicts a faint 
aurora (1100 rayleighs). However, in addition to 
explaining the intensity level, a successful theory 
also has to account for the time variations of the 
intensity. On a purely static picture, it is not 
possible to understand why the intensity maximum 
does not come in the beginning of the night. Since 
high wind speeds of the emitting cells have been 
deduced the dynamics of the airglow emission into 
the photochemical theory will now be introduced to 
explore the possibility of nocturnal time variations 
as observed. 

3. Meterological Approach 

As long as airglow is considered as a steady, static, 
more or less homogeneous emission from a layer at 
about 100 km, it is natural to try to explain it in 
terms of a pure photochemical process solely gov- 
erned by the kinetic properties of the atmosphere. 
With the introduction of a dynamical picture in 
which the airglow emission is concentrated in cells 
that move relative to the earth's surface, one is led 
to ask whether the very dynamics of the emissions 
have anything to do with the excitation of the 
phenomenon. 

At least two avenues lead from here. First, one 
may ask whether the dynamics of the upper atmos- 
phere can modify the emissions caused by the photo- 
chemical reactions to a sufficient degree to account 
for the observed morphology of the airglow. This 
may be termed the meteorological approach. Second, 
one may go a step further and inquire into the possi- 
bility of using the dynamics also for the excitation of 
the oxygen atoms. In this case, the energy of the 
environmental effects is used as in the next section. 

In light of the increased interest of recent years in 
the study of solar-weather relationships and in the 
possibility of a coupling between the troposphere and 
high atmosphere (including the lower part of the 
ionosphere), it seems appropriate to consider in some 
detail the meteorological approach to the airglow 
problem. 

The existence of cells in the airglow emission is 
probable, and from a meteorological point of view, 
they may be considered as "clouds." The very 
concept of clouds is synonomous with instability, and 
the first question is whether one would expect insta- 
bility in the high atmosphere similar to that in the 
troposphere. Fleagle [14] has recently shown that 
disturbances tend to grow rapidly in the mesodecline, 



4 If the brightness, B, is measured in units of 10 6 quanta/cm2 second steradian, 
then the brightness in rayleighs (ff) is 4irB. 



at least in polar regions. Also, it should be remem- 
bered that the outstanding feature of the troposphere 
in this connection is the lapse-rate — the decreasing 
temperature with height and the resulting small 
Richardson number (Richardson [15]) which is I he 
criterion that turbulence in a gravitational field will 
increase. 5 (In contrast, the high Richardson number 
in the stratosphere points to stable conditions.) The 
Richardson number is given as: 






(9) 



where V is the adiabatic lapse rate (=—), 

c p = specific heat at constant pressure, 
g= acceleration of gravity, 
T= temperature, 

##= horizontal velocity component of the wind, and 
z is in the vertical direction. 

As an example, the Richardson number for the 
lower troposphere, assuming the following character- 
istics: 

c p = 1.003 j/g°K, 

T=270°K, 



dT 
dz 

dvn^ 
dz 



= -5X10- 7O K/cm 



:1.9X10- 2 /sec, 



is 



R t (troposphere) = 0.98. 



The only other region in the atmosphere where a 
decreasing temperature with height is found is in 
the mesodecline, at the top of which the airglow 
phenomenon is observed. Data, are not so easily 
available to compute the Richardson number there, 
but the vertical change in the wind velocity is of 
the same order of magnitude as in the troposphere 
and instability would be expected. 5 In other words, 
convection is likely to occur and it should reach 
roughly to the region of minimum temperature, the 
mesopause, and this is roughly where we find the 
green airglow. 

If the emission changes can be correlated with the 
convective motions, then, on this picture, the airglow 
cells are to be considered as the top part of the 
convection elements in the mesodecline. It is a 
priori not possible to predict in any detail how the 
penetration of air masses into the atomic oxygen 
concentration peak will affect the green emission, 
and only the conjecture is made here that this is the 
reason for the time variation of the airglow inten- 
sities. 



a It is generally agreed (Godske, Bergeron, Bjerknes, and Bundgaard [16]) 
that a Richardson number somewhat less than unity points to instability. But 
the exact critical value of the Richardson number is not known. 
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Now, the mesopause is a rather shallow minimum. 
One even thinks of it as an isothermal region situ- 
ated at around 85 km. This is below the airglow 
height and this looks, at first sight, somewhat incon- 
sistent. Two points should, however, here be taken 
into consideration. First, as Nicolet [10] points out, 
the height of the mesopause is not known to the 
nearest kilometer. Nicolet gives as the best esti- 
mate 85 ± 10 km. Second, if the region around 90 
km is nearly isothermal, one would expect the con- 
vection elements to penetrate up through this layer, 
similarly as one observes eddies in the region just 
above the tropopause. 

The interpretation of our data as wind implies that 
we have a component of ordered motion, a sort of 
"jet-stream," near the mesopause. On the other 
hand, this motion is certainly superposed on turbu- 
lence. The Reynold's number is 



R- 



Lv 



where I is a characteristic length, v is a character- 
istic velocity, and v is the kinematic viscosity. With 
U= 10 4 cm/sec and p=4X10 4 cm 2 /sec, a Reynold's 
number R e =0.25 L (L in cm) is found. It is not 
obvious which characteristic length to put in, but 
if the boundary layers of the cells involved in the 
turbulence have a thickness as small as 1 km, 
R e =2.5X10 4 is found. Convection will tend to 
inhibit the turbulence, but since this is a region to 
which the convection only just penetrates, the char- 
acteristics of turbulence should probably manifest 
themselves. Evidence of this may be the rapid 
changes in luminosity often observed in the airglow 
cells. 

On this picture, the maximum airglow intensity 
during a night may now come at any time, depending 
on the penetration of the instabilities into the maxi- 
mum layer of the atomic oxygen distribution. 
Since the intensity will vary as [n(0)] 3 rather great 
variations may be expected, and due to the nature 
of turbulence relatively erratic time variations are 
not excluded. The green airglow is hence to be con- 
sidered as a manifestation of "the weather in the 
mesodecline." 

Obviously, too little is known about the dynamics 
of the regions involved to draw any definite conclu- 
sions, but this approach may be worth further study 
as more data on this part of the atmosphere become 
available. 

Attention is called to a study of 5577 by Tohmatsu 
[17] in which the Chapman reaction is assumed and 
the localized variations in atomic oxygen concen- 
tration are invoked to account for the observed 
patchiness of the night airglow, as well as the large 
scale temporal and geographical charges. 

4. Environmental Effects 

The fact that the 5577 airglow emission occurs in 
the ionosphere will now be considered. The electron 
density is small compared to the neutral gas density, 



but even in such a slightly ionized plasma new effects 
are likely to occur. This is the more probable since 
the plasma is pervaded by and moves in the earth's 
magnetic field. 

The ionosphere has a general similarity to a gas 
discharge tube, and it is reasonable to inquire 
whether the conditions that lead to discharges in the 
laboratory are present in the ionosphere and will 
lead to luminous discharges that can be associated 
with airglow. 6 It is characteristic of gas tube dis- 
charges that they can be explained only on the 
assumption that the electrons present are able to 
absorb energy selectively from an electric field 
between elastic collisions with heavy particles in the 
plasma. The essential feature is hence the presence 
of an electric field capable of accelerating electrons 
until they have enough energy to excite oxygen 
atoms. This leads to conditions in the plasma 
radically different from equilibrium and charac- 
terized by an electron gas temperature greatly in 
excess of the heavy particle temperature. 

Several attempts have been made to formulate the 
problem. Chamberlain [18] postulated the existence 
of an electric field and determined the energy dis- 
tribution curves and drift velocities of the electrons. 
He was able to show that auroral rays could be 
explained on such a discharge theory. Alfven [19] 
and Wulf [20, 21] also discussed the mechanism by 
which the electric field could be set up. The problem 
includes two main facets: (a) The existence of an 
electric field capable of accelerating electrons, and 
(b) the heating of the electron gas and the excitation 
of the oxygen atoms. 

(a) Production oj an electric field. 

Since it is implied here that the dynamics of the 
100-km region are of prime importance to an under- 
standing of the green emission, those attempts to 
establish an electric field which involve extrater- 
restrial streams (see Alfven [19]) will not be con- 
sidered. This stage of the discussion is more in line 
with Wulf [21], inasmuch as the velocity V c of the 
airglow cells as they move across lines of force of the 
earth's magnetic field B will be used. In doing so, 
an electric field is induced 



E lnd = V c XB (emu) 



(10) 



This field will tend to produce an electrostatic 
separation of charges until the polarization field 
Epoi thus created exactly balances Ei nd . Hence, in 
the frame of reference of the cell, there is no net 
electric field. Or, expressed in another way, the 
polarization field gives the charged particles (irre- 
spective of sign) a drift velocity Y d in the direction 
of the motion V c , such that V d =V c , The plasma 
therefore moves "downstream" as an entity, and 



EXB 



B 2 



(11) 



6 It is a pleasure to acknowledge the many stimulating discussions on this 
topic with Dr. E. N. Thomas who first drew our attention to this similarity. 
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In the absence of collisions, the electrons do not 
sec any field and a heating of the electron gas is not 
possible. When t he collisions are taken into account, 
I he problem becomes quite difficult, and it is not at 
all obvious I lint the electrons can draw energy from 
the Held. Wulf [21] avoided this basic difficulty by 
letting the circuit be closed through the F 2 layer 
where the discharge took place. 7 This would be 
physically understandable, but it is known the air- 
glow emission comes where the wind speed V c is 
measured, and this rules out the original Wulf 
approach. It may be recalled that Wulf used eq 
(10) and found that with a value of 0.5 gauss for the 
earth's magnetic field and a wind speed of 50 m/sec 
a gradient of 

#=0.5X5X10 3 = 2.5X10" emu 

or £~10~ 5 volts/cm 

would he created. As mentioned, this potential drop 
is not available for electron acceleration in this case, 
since an external path over which to discharge is 
needed. 

Another possibility remains. Roach et al. [5] 
found evidence of a rotational velocity of the order 
of 400 m/sec in the airglow cells. If this can lead to 
polarization effects in the cells, stronger fields (per- 
haps 10" 4 to 10" 3 v/cm) might be expected to be set 
up between adjacent cells due to the great differen- 
tial peripheral speeds. But the 4 whole question is for 
the time being wide open, and no conclusion should 
be drawn before this problem lias been thoroughly 
analyzed. The best for the moment is to assume the 
existence of the field and see what the results might 
be. 

(b) Acceleration of electrons and excitation of oxygen. 

If the same point of view as Chamberlain's [18] is 
taken, the question may be asked how a postulated 
electric field will affect the electron energy distribu- 
tion. According to Chamberlain a gradient of order 
of magnitude between 10~ 5 and 10 _1 v/cm is sufficient 
to account for the acceleration of electrons to an 
energy capable of exciting the oxygen atoms in 
auroral rays. Now, it should be remembered that 
this conclusion is reached on the basis of deriving the 
electron energy distribution by a simplified theory 
(Smit [22]). The complete treatment of the problem 
from the appropriate Boltzmann equation remains 
to be tackled. 

The intensity of the 5577 emission is, according to 
Chamberlain, given by 



Q(5577) = 



A, 



A 32 +^-31 



n(OI)n e S v 



io 6 



where collisional deactivation lias been 
and Q is in rayleighs. 



(12) 
'gleeted. 



7 At that time, it was believed that the green airglow originated al a height of 
approximately 200 km. 



Using a thickness of the emission layer equal to a 
scale height (27=6.5 km), the values given in the 
right side of figure 2 (n e =10 2 to 10 5 ) are obtained; 
and with // = 100m(to compare with Chamberlain's 
auroral conditions) the values in the left side of 



figure 



(for n e =l0 4 to 10 s ). An energy e 



/'A 



where A is the mean free path) of approximately 
10" 1 v is necessary to explain the airglow intensities. 
The extremely sharp dependence of the emission on 
this parameter (e) may have an important bearing on 
the intensity variations in time and space of the 
airglow pattern, since 4 inhomogeneities in E are likely 
to occur. 




10' 10" 

6 = E X (volt) 
Figure 2. Predicted intensity of [OI] 5577 from eq {12), 



The airglow is a normal feature of the "quiet" 
atmosphere in the 100-km region, and it would here 
seem reasonable to compare it with other quiet- 
condition phenomena of electromagnetic nature. 
The daily variations in the geomagnetic field are such 
a phenomenon. The variations can be explained in 
terms of a current system in the upper atmosphere, 
and rocket observations (Cahill [23]) show the exist- 
ence of currents at 95- to 110-km altitude. The 
question may therefore be asked whether electrons 
could be sufficiently accelerated in the electric field 
associated with this current system. On the night 
side of the earth where the airglow is observed, the 
current between equator and 40° north is approxi- 
mately 40,000 amp (Chapman and I hut els [24]). 
The width of the current sheet is hence 40X110 
km -4,000 km. The thickness / of the sheet is not 
too well known so the current density j is written as 



3=- 



M :,(>!)( I r.'.t 
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where the area ^=ZX4X10 8 om 2 (I in cm). Using 
Ohm's law, 

j = aE (13) 

the corresponding electric field E is found provided 
the anisotropic electrical conductivity a is known. 
On magnetically disturbed days additional intense 
current systems at high latitudes can be estimated 
Matsushita [25]). Under quiet conditions (Mat- 
suchita [26]) the electric field is estimated to be 
approximately 3X10 3 emu = 3X10~ 5 v/cm. This 
is on the low side of what figure 2 indicates is nec- 
essary to produce the observed emission. On the 
other hand, the value 3X10 -5 v/cm is clearly a 
lower limit, since it is deduced using quiet geomag- 
netic conditions. The existence of magnetic dis- 
turbances means a greater electric field, and on the 
discharge theory this may be sufficient for the neces- 
sary heating of the electron gas. 

It is hardly necessary to stress that this per se is no 
theory for the airglow, but the calculations indicate 
that it may be worth looking for an explanation of 
the X5577 emission along these lines. 

5. Summary and Conclusions 

The preceding paragraphs have examined existing 
theories and hypotheses of the green oxygen emis- 
sion and presented promising alternatives. The 
sparsity of data available for some of the coefficients 
and parameters involved makes numerical calcula- 
tions difficult, but some general conclusions can be 
drawn. 

It is believed that "a modified photochemical 
theory" is capable of explaining the observed 5577 
emission. Tin* complex dynamics of the 90- to 100- 
km region are not well known, but from what can 
reasonably be deduced, the "weather conditions" in 
the mesodecline are compatible with the observed 
variations of the 5577 emission in space and time. 

Similarly, it is possible that the emission can be 
understood on the basis of "the environmental ap- 
proach." In particular, the steepness of the curves 
in figure 2 indicates that small changes in electron 
density and/or the electric field (which probably occur 
in the low flayer at night) will lead to rapid changes 
in brightness, and may explain space and time varia- 
tions of the emission. 

From the present study, one mechanism cannot 
be pointed to as the only possibility. This is per- 
haps not a very happy situation, but it may be that 
Nature is using two different mechanisms to generate 
the observed emissions. 



The research reported in this paper has in part been 
sponsored by the Geophysics Research Directorate 
of the Air Force Cambridge Research Center, Air 
Research and Development Command, under Con- 
tract AF 19(604)-2140 to the High Altitude Observa- 
tory. 
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